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I critically review the status of computations of threshold pion production in 
C J , proton— proton collisions in the framework of effective field theory approaches or 

variants thereof. I also present the results of a novel diagrammatic scheme. 



1 The problem 



> . 

*n i Over the last years, very precise data on pion, rj and rf meson production in 

proton-proton collisions in the threshold region have been obtaiaed at IUCF, 
CELSIUS and COSY (for a review, see the talk by H.O. MeyerH). The basic 
process is pp ~ > pNM, where N denotes the nucleon and M a pseudoscalar 
meson, in our case the ir°, tt + , r\ or the rf . At the respective threshold, the 
produced meson is soft, i.e. has vanishing three momentum. Consequently, in 
the case of the pions, which are believed to be the Goldstone bosons related to 
the spontaneous chiral symmetry breaking QCD is assumed to undergo, this 
process appears to be a good testing ground for chiral perturbation theory 
^ \ methods. To be specific, at threshold one has only S-waves and thus the 

pertinent T-matrix for the specific case of ir° production is parametrized in 
terms of one single amplitude, 

X' T^(pp^ppTt°) =A(iai-ia 2 + ai x ff 2 )'P ■ (1) 

The (71,2 are the spin-matrices of the two protons. The amplitude is a pseu- 
doscalar quantity and due to the Pauli principle, we are dealing with a 3 P — > 
1 Sq s transition (where the V refers to the pion angular momentum). The 
value of the proton cm momentum to produce a neutral pion at rest is given 

by 

\p\ = a/M^ (m + Mjl) ~ s/~M^m = 362.2 MeV , (2) 

with m = 938.27 MeV the proton and M n = 134.97 MeV the neutral pion 
mass, respectively. Obviously, \p\ vanishes in the chiral limit of zero pion 
mass. Therefore the soft-pion theorem which requires a vanishing threshold 
T-matrix in the chiral limit M n — is trivially fulfilled (as long as A does not 
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become singular). Stated differently, there is no low-energy theorem for the 
threshold amplitude A. Note also that the value of the threshold momentum is 
already relatively large. Furthermore, it is known that in the threshold region 
the strong final-state interactions (FSI) govern the energy dependences In 
fact, in the approaches I will discuss in what follows, the complete amplitude 
is written as 

rp rplSl rpProd rpFSl 

where ISI denotes the initial-state interaction and the microscopic approaches 
are applied to the production amplitude T Prod . Clearly, such a separation 
induces a priori some model dependence, as I will discuss later. Let me first 
describe the existing EFT approaches to calculate/constrain the production 
amplitude. 



2 Heavy fermion techniques at work 

Most work so far has focused on the production of neutral pions simply because 
the production amplitude involves, besides many other processes, isoscalar 
pion-nucleon scattering (see fig. lb), which is known to be suppressed due to 
chiral symmetry. However, it should be made clear frorn_the start that there 
must also be important short distance physics, see e.g. a and as depicted in 
fig.lc. The challenge of working out the pion exchange terms in chiral pertur- 
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Figure 1: Different mechanism contributing to the meson production operator, a), b) and 
c) are referred to as the direct, rescattering and heavy meson— exchange terms, in order. 

bation theory has been been taken up by two groups, one at Seattle^ and the 
other at South CarolinaB Their calculational schemes are similar. To be defi- 
nite, I concentrate in what follows on the work of the Seattle group. They use 
the effective pion-nucleon Lagrangian to second order in small momenta (i.e. 
tree level graphs), the appearing low-energy constants (LECs) have already 
been determined from pion-nucleon scattering datau The nucleons are treated 
as very heavy spin-1/2 fields with four-momentum 

Pn = mv^ + fc M , (4) 
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with Vfj, the four-velocity subject to the constraint v 2 — 1 and v ■ k -C m. To 
leading order, the nucleon propagator is S(k) = i/(v ■ k) and one can set up a 
consistent power counting scheme in terms of small momenta, the pion mass 
and inverse powers of the nucleon mass. However, their calculations involve a 
modified counting scheme due to the abovementioned fact that the momentum 
does not scale like the pion mass, so the heavy fermion propagator for a small 
residual momentum fc M includes the first kinetic energy correction 

(5) 



v-k k + k 2 /2m 

since in the rest frame = (1,0). The necessity of such a modification follows 
from the fact that the heavy baryon formalism can not cope with external 
momenta as large as \p\ ~ \J mM v . Let me demonstrate the problems of 
the heavy baryon approach for a simple (but generic) example. Consider a 
Feynman diagram which involves a propagating nucleon after emission of the 
real tt°. I show that this nucleon propagator can not be expanded in powers 
of 1/m in the usual way. Let u M = (1,0) be the four-vector which selects the 
center-of-mass frame. The four-vector of the propagating nucleon is m v + k 
with fc M = (—M n /2,p) and k 2 = —mM^. Starting on the left hand side with 
the correct result based on a fully relativistic calculation and then performing 
the usual 1/m expansion of the heavy baryon formalism, one has 

- j_ = i = — f r = -— f ( _ ir (6) 

M n v-k + k 2 /2m vk^\2mvkj M„ ^ y ' w 

71 1 n=0 n=0 

One sees that infinitely many terms of the 1/m-expansion contribute to the 
same order. The resulting series does not even converge and oscillates between 
zero and twice the correct answer. The source of this problem is the extreme 
kinematics of the reaction NN — > NNir with \p\ ~ \/mM„. In that case the 
leading order operator = iv ■ d and the next -to-leading order operator 
0( 2 ) = —d ■ d/2m lead to the same result, here, iM^/2. This problem is 
merely related to "trivial" kinematics. Therefore, one either has to modify the 
propagator as described above or calculate fully relativistically, as discussed 
in the next section. Space forbids to discuss in big detail the results of these 
studies, I rather concentrate on the most intriguing ones. First, one finds 
that the rescattering contribution (fig. lb) comes out with an opposite sign to 
what one gets in meson-exchange models, i.e. it interferes constructively with 
the direct production in the Jiilich meson-exchange model u and destructively 
in the chiral framework, respectively. Note that while there is still debate 
about the actual numerical treatment and the ensuing size of the rescattering 
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contribution in the chiral perturbation theory approaches^ the sign difference 
with the meson-exchange-model can be considered a genuine feature. This 
point was addressed in rem. It was argued that the treatment underlying the 
isoscalar pion-nucleon scattering amplitude and the related transition operator 
for the process NN — > NNtt in the chiral framework is not yet sufficiently 
accurate and thus the resulting rescattering contribution should be considered 
an artifact of this approximation. Clearly, this does not mean that chiral 
perturbation theory is invalid but rather that higher order (one loop) effects 
need to be accounted for. Second, to account for the data, one has to include 
short distance physics, like e.g. heavy meson exchanges (see fig.lc) and also 
meson-exchange currents, in particular the anomalous irpuj vertex. A typical 
result is shown in fig. 2, which also shows that there is still a sizeable uncertainty 
related to e.g. the treatment of the isoscalar rescattering. 
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Figure 2: Cross-section for pp — > ppiz as function of the pion momentum r) in units of M w 
for two NN potentials (Argonne V18 and Reid93) and two parameter sets for the isoscalar 
7riV amplitude (ste and cl). Figure courtesy of Bira van Kolck. 

The first one-loopxalculation (to third order in small momenta) was per- 
fomedjpy Gedalin et alE3. They work in the conventional heavy baryon frame- 
work EH using the lowest order fermion propagator (see the left-hand side of 
eq.(|^)). Some typical one-loop graphs are shown in fig. 3. There are, of course, 
many more diagrams, most of them giving rise to mass and coupling constant 
renormalization. At this order, there appears also a four-nucleon-pion contact 
term with a LEC, called d\. This LEC could in principle be fitted from one 
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Figure 3: Class of one loop graphs involving the -kit interaction. Dashed/solid lines denote 
pions / nucleons . 



total cross section point. The authors of ref0 estimate its values from reso- 
nance saturation, specifically by p and w exchange. This procedure induces, 
of course, some model-dependence. (The status of resonance saturation in the 
one-nucleon sector is discussed in refcl) Therefore, the total cross section can 
be predicted without any free parameters and it agrees quite nicely with the 
data in the threshold region. However, there are a couple of loopholes with this 
result. First, as discussed before, using the leading order nucleon propagator, 
one can not expect a convergent series. This is reflected in the results of reftJ, 

^irescattcring . ^loop . ^direct ~ 5 • 2 ' 1 (7) 

according to figs. lb, 3, la, respectively. Without the loop contribution, the total 
cross section is underestimated by a factor of about 2.5. One can pin down 
this problem of the bad convergence more specifically. The class of loop graphs 
shown in fig. 3 is directly proportional to the scalar form factor of the nucleon 
at t = M^m ~ 0.1 GeV 2 . In the heavy fermion approach, it leads to a very 
large enhancement of the production amplitude. If one, however, calculates 
these loop diagrams fully relativistically and evaluates the scalar form factor, 
its contribution is very smallt 2 ] This simply reflects the fact that a one-loop 
calculation at third order is not sufficiently accurate, even the shift of the scalar 
form factor to the Cheng-Dashen point at t = 4Af 2 ~ 0.02 GeV 2 is off by a 
factor of two. Stated differently, the relativistic calculation shows that higher 
order 1/m terms are large. Similar remarks apply to the full isoscalar pion- 
nucleon scattering amplitude. From these results I conclude that heavy baryon 
chiral perturbation theory is presumably not the appropriate framework to gain 
a deeper understanding of tharole of chiral dynamics in the process pp — > ppir . 
As already pointed out in rem, a more detailed study of the reaction pp — ► dir + 
could, however, pave the way to a deeper understanding of the pion dominated 
part of the transition operator (simply because short-distance physics plays 
a much smaller role in this process). Finally, I remark that the treatment 
discussed so far involves an unavoidable ambiguity due to the exchanged pion 
(say in the rescattering graph, fig. lb) being off-shell. For a given prescription 
of defining the pion field, one can perform an off-shell extrapolation but this 



will always be model-dependent - in quantum field theory only on-shell matrix 
elements and transition currents can be calculated. 

3 Diagrammatic approach 

Based on the observation that the heavy baryon framework is not expected to 
converge, in refO a different approach was pursued. Consider first the process 
pp — > ppir . Approximating the near threshold T-matrix by the T-matrix 
exactly at threshold one gets for the unpolarized total cross section 

<7tot(T lab ) = \A\ 2 J dWKF(W)F p [W) , (8) 

where the flux and three-body phase space factors, denoted as KF(W), can be 
approximated by an analytical expression which is accurate within a few per- 
cent in the threshold region. F p (W) is the correction factor due the final-state 
interaction, with W the final-state invariant di-proton mass. It is evaluated 
in the effective range approximation, 

F P (W) = {l + ^(a p + r p )Pl + ^fP^ 1 , 0) 

with P 2 = W 2 — 4m 2 and a p (r p ) the pp scattering length (effective range) 
including electromagnetic corrections. This is of course a very strong assump- 
tion but it allows one to explain the energy dependence of the experimental 
total cross sections very accurately in terms of a single constant amplitude 
A. At the end of this section, I derive this particular treatment of the FSI 
from an effective field theory (EFT) approach. Separating off the final-state 
interaction in that way, one can then pursue a diagrammatic approach to the 
(on-shell) production amplitude A. This allows one to investigate in a sim- 
ple fashion the role of one-pion exchange and chiral loop effects together with 
shorter range exchanges due to heavier mesons. In a similar fashion, one can 
investigate the other processes pp — > pmr + , pp — > pprj and pp — > pprj 1 (see 
below). To be specific, consider the ppir reaction. Assuming the form eq.(|J), 
the S-wave threshold amplitude can be extracted from the data, 

A {cmp] = (2.7 -* 0.3) fm 4 , (10) 

as shown by the solid line in fig. 4. The imaginary part is due to the 3 Pq pp phase 
shift taken at the threshold energy in the lab frame, T^ h = 279.65 MeV, where 
5( 3 Pq) = —6.3°. Thus the imaginary part Im„4 is about —1/9 of the real part 
Re A and contributes negligibly to the total cross section near threshold. This 
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Figure 4: Fit to the total cross section for pp — ► ppn° as described in the text (solid line). 
The data are from IUCF (boxes) and CELSIUS (crosses). The dashed line is result of the 
diagrammatic approach. 



number can be well understood in terms of chiral ir° exchange (including chiral 
7T° rescattering) and heavy meson (u>, p° , rf) exchanges based on a relativistic 
Feynman diagram calculation with all parameters being fixed from reliable 
methods, like forward NN dispersion relations for the u and a dispersion- 
theoretical analysis of NN — > irir for the p. Interestingly, in the relativistic 
approach the rescattering contribution interferes constructively with the direct 
production term (cf. fig. la), in contrast with the heavy baryon approach and 
in agreement with the meson-exchange models. One can also evaluate some 
classes of one-loop graphs and finds that they lead to small corrections of 
the order of a few percent. Therefore chiral loops do not seem to play any 
significant role in the processes NN — > NNtt, which are dominated by one- 
pion exchange and short-range physics. I remark also that both the long range 
7T° exchange and the short range vector meson exchange lead to contributions 
to the threshold amplitude A which do not vanish in the chiral limit M t — > 0. 
There is no chiral suppression of the reaction pp — ► ppir compared to other 
NNtt channels. In all cases the respective threshold amplitudes are non-zero 
(and finite) in the chiral limit. This is in contrast to the widespread believe 
that pp — > ppir is suppressed for reasons of chiral symmetry. 

Within the same approach, one can investigate the threshold behavior of 
the process pp — > pnir + . It is given in terms of A and the triplet threshold 
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amplitude B, 



It? 1 (PP P^tt + ) = (tax - i a 2 + 3\ x ct 2 ) • p - \/2Si(CTi + ct 2 ) • p ■ (11) 
v2 

Here, the second term refers to the 3 Pi ^ 3 Sis transition. Note also that the 
singlet transition is suppressed because of the large singlet scattering length. 
From the IUCF data, one can determine the empirical value of the triplet 
amplitude, B cmp — (2.8 — il.5) fm 4 . The large imaginary part is related to the 
strong ISI in the 3 Pi entrance channel, which naturally can not be explained 
by tree graphs only. The value of the 3 Pi phase at the threshold energy for 
pp — > pmr + is —28.1°. The corresponding real part Re B is well reproduced 
by chiral one-pion exchange and short-range vector meson physics. For the 
same parameters as used in the study of pptr , one gets Re B — 2.74 fm 4 . A 
more detailed discussion of this chaxmel and the problems related to the large 
imaginary part can be found in reflij. 

Let me now turn to r\ production. The threshold matrix element takes the 
form 

Tth (PP -» PPV) =C{ia 1 ~ia 2 +S 1 xa 2 )- p. (12) 

with C the (complex) threshold amplitude for ^-production. The ^-production 
threshold is reached at a proton laboratory kinetic energy TJ*^ = M n (2 + 
M v /2m) = 1254.6 MeV, where M v = 547.45 MeV denotes the eta-meson 
mass. In the case of ^-production near threshold it is also important to 
take into account the rjp final-state interaction, since the ryiV-system inter- 
acts _cather strongly near threshold. In fact a recent coupled-channel anal- 
ysis E3 of the (7riV, ?yiV)-system finds for the real part of the r]N scattering 
length Kea V N — (0.717 ± 0.030) fm. For comparison, this value is about six 
times larger than the ir~p scattering length, a^- p — 0.125 fm, as measured 
e.g. in pionic hydrogen. In refO, it is assumed that the correction due to 
the S-wave rjp FSI near threshold can be treated in effective range approxi- 
mation analogous to the S-wave pp FSI. The further assumption is made that 
the FSI in the pp subsystem and in the two rjp subsystems do not influence 
each other and that they factorize. The corresponding form of the unpolar- 
ized total cross section in terms of the S-wave amplitude C the various FSI 
functions F p (W), F rj (s n ) can be found in refEil. Note that the r\p FSI function 
is complex-valued (because the rjN scattering length has a sizeable imaginary 
part). From the CELSIUS data, one finds for the modulus of the threshold 
amplitude \C\ = 1.32 fm 4 . The resulting energy dependent cross section from 
threshold up to Xi a b = 1375 MeV is shown in fig. || together with the data from 
CELSIUS. It is rather astonishing that one can describe the total cross section 
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Figure 5: The eta-production cross section crtot(pp — > ppv) as a function of Ti a t,. The data 
are taken from CELSIUS. 



data up to 100 MeV above threshold with a constant threshold amplitude C 
and a simple factorization ansatz for the three-body FSI. The relativistic Feyn- 
man graphs contributing to the threshold amplitude C can be easily evaluated. 
To account for the strong r]N S-wave rescattering, which is often attributed to 
the nucleon resonance Sii(1535), one has to introduce a local NNrjrj contact 
term, C^n = K N (x) N (x)t] 2 (x) . Its strength K can be obtained from fitting 
the scattering length a V N- With g V N = 5.3, which is close to flavor SU(3) 
estimates and determinations from boson-exchange models, one can exactly 
reproduce the value of C. It is worth pointing out that p° exchange is the 
dominant contribution, because it is enhanced by factors of M v /m ~ 0.6 quite 
in contrast to the neutral pion case, where ui and 7r exchange are the dominant 
mechanisms. One can also predict the threshold S-wave amplitude T> for the 
process pn —> pnrj. One finds \T>\ — 1.8 fm , somewhat below the experimen- 
tal value, |2?| cmp = 2.3 fm 4 , as determined from the recent data on quasifree 
production off the deuteron measured at CELSIUS. 

At COSY, rj production very close to threshold has been measured for the 
first time, at much smaller excess energies than the few data points from the 
now defunct SATURNE facility. These data seem to follow three-body phase 
space with little indication of FSI. Note, however, that so close to threshold 
the energy determination is difficult and thus there is a sizeable uncertainty 
in the value of the excess energy (for more details, see refo). Within the 
diagrammatic approach, the empirical S-wave amplitude can be obtained with 



9 V 'NN (1 - 1.28e) = 1.12 , 
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(13) 



with g v i jv at the so far undetermined rf — TV coupling constant and e the pseu- 
doscalar to pseudovector mixing parameter. Since the rf is not a Goldstonc 
boson, one is not forced to use only derivative couplings. Interestingly, only 
the tensor interaction of the p-exchange (~ k p ) is sensitive to the parameter 
e. Combining this result with constraints from deep inelastic lepton-nucleon 
scattering!!!] (the spin content of the nucleon) , one infers 

g v , NN = 2.5 ±0.7, £ = 0.4 ±0.1. (14) 

It remains to be seen whether other ^'-production processes (e.g. photopro- 
duction 7p — > rfp) are consistent with these values. 

Finally, I want to give an elementary derivation of the final-state interac- 
tion correction factor F p (W) = [1 + a p P 2 ]~ 1 , P 2 = W 2 /4:~ m 2 , in the scatter- 
ing length approximation (i.e. for the effective range r p = 0). Close to thresh- 
old all final state three-momenta are small and therefore one can approximate 
both the meson production process TV TV — ► NNir and elastic TV TV — * TV TV 
scattering by momentum independent contact vertices proportional to A and 
the scattering length a p , respectively. Consider first low energy TV TV-scattering 
in this approximation. The bubble diagrams with 0,1,2,. . . rescatterings can be 
easily summed up in the form of a geometric series, 



v 



- I { i ) 3 pi + i 0+ _ P +- = 'h + i^ + - = ' (15) 



using dimensional regularization (in the MS'-scheme) to evaluate the (vanish- 
ing) real part of the loop integral. Obviously, the sum of these infinitely many 
loop diagrams is just the unitarized scattering length approximation which 
leads to tanSo(W) = a p P» . Next, consider in the same approximation me- 
son production followed by an arbitrary number of iVTV-rescatterings in the 
final-state. Again, these loop diagrams can be summed up to, 



A 



1 2 CLp -tj(s 



and taking the absolute square, 

A 



1 — i a p P* 



\A\'< 



1 ± a 2 p P? 



(16) 



(17) 



one encounters the final-state interaction correction factor F p (W) = [1 + 
a p (W 2 /4 — to 2 )]" 1 in scattering length approximation. Since the scattering 
length is much bigger than the effective range parameter for iV TV-scattering 
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(a p >> r p ) one has already derived the dominant effect due to the FSI. Of 
course, in order to be more accurate one should eventually go beyond mo- 
mentum independent contact vertices. The main point, however, I want to 
make here is that the FSI correction factor F V (W) (for r p = 0) has a sound 
foundation in effective field theory. 

Of course, this approach also has some drawbacks. First, in view of the 
recent polarization measurements at IUCF and RCNP, one should go beyond 
the S-wave approximation. Second, it is not obvious how to systematically 
improve the calculations and third, the treatment of the FSI might be too 
simplistic (for a recent discussion, see e.g. reffi3). 

4 Outlook 

I given a short review of EFT or EFT-inspired approaches to pion (as well as 
r\ and ?/) production in proton-proton collisions. The much discussed reac- 
tion pp — > ppir does not seem to offer a testing ground for chiral dynamics 
due to large uncertanties related to chiral pion-exchange, the short distance 
physics and the strong ISI/FSI. Even the novel EFT scheme due to Kaplan, 
Savage and Wise EH is (in its present formulation) not able to cope with pion 
momenta as large as in this reaction. On the other hand, there exist also the 
rather successfull meson-exchange models based on effective Lagrangians from 
the Jiilich (and RCNP) group(s). While these work fairly well, there is no 
systematics means of improving them and also, they employ highly model- 
dependent meson-nucleon form factors, which have no sound physical basis. 
Since a consistent EFT approach to deal with these processes does not seem 
to be on the horizon, it might be worthwhile to try to constrain the meson- 
exchange models with EFT ideas and eventually gain a better understanding 
of the meson-nucleon interaction regions. Also, more effort should be invested 
in studying the process pp — > dir + since it is much less contaminated by short- 
distance physics. Theory still has a long way to go to catch up with the 
tremendous precision obtained in the recent experiments. 
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